phase transition calculation